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We calculate the ion and composition distributions around colloid particles in an aqueous mix- 
ture, accounting for preferential adsorption, electrostatic interaction, selective solvation among ions 
and polar molecules, and composition-dependent ionization. On the colloid surface, we predict a 
precipitation transition induced by strong preference of hydrophilic ions to water and a prewetting 
transition between weak and strong adsorption and ionization. These transition lines extend far 
from the solvent coexistence curve in the plane of the interaction parameter \ (or the temperature) 
and the average solvent composition. The colloid interaction is drastically altered by these phase 
transitions on the surface. In particular, the interaction is much amplified upon bridging of wetting 
layers formed above the precipitation line. Such wetting layers can either completely or partially 
cover the colloid surface depending on the average solvent composition. 
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I. INTRODUCTION 

Extensive efforts have been made to understand the 
interaction among ionized colloid particles in a solvent 
P-Q , because they form model crystal and glass at high 
densities. Recently, considerable attention has also been 
paid to the effect of preferential adsorption of one of the 
components in a mixture solvent. Several groups [rjl-floj 
observed aggregation of colloidal particles near the co- 
existence curve in one-phase states of a binary mixture 
of 2,6-lutidine and water. For small ionization the ad- 
sorption of lutidine was increased in water-rich states, 
while for large ionization that of water was increased in 
lutidine-rich states. This means that the colloid surface 
can either repel or attract water, depending on the degree 
of ionization. It is worth noting that polyelectrolytes are 
often hydrophobic without ionization, but becomes effec- 
tively hydrophilic even at low ionization (ll| - [l3j . 

At high colloid concentrations, a flocculated phase rich 
in colloids emerges [6], which changes from gas, liquid, fee 
crystal, and glass with increasing the colloid concentra- 
tion Q. Such aggregation was claimed to be a result 
of a true phase separation in ternary mixtures Q. A 
microscopic theory by Hopkins et al. [l4j indicated the 
aggregation mechanism. They treated neutral colloids 
coated by a thick adsorption layer rich in the preferred 
component in one-phase environments rich in the other 
component. They found that this adsorption is intensi- 
fied near the coexistence curve, strongly influencing the 
colloid interaction. 

In understanding these experiments, however, atten- 
tion has not been paid to the selective solvation (hydra- 
tion for aqueous mixtures) among charged particles (ions 
and ionized parts on the colloid surface) and polar sol- 
vent molecules [15]. The solvation effects have not yet 
been adequately investigated in soft materials, but they 
can influence the phase separation behavior profoundly 
and even give rise to a new phase transition |16H19| . The 
solvation chemical potential /i so i(</>) of a hydrophilic ion 
stems from the ion-dipole interaction and strongly de- 



pends on the ambient composition of water. In electro- 
chemistry, the following chemical potential difference has 
been measured 120, EH : 



(1.1) 



between two coexisting phases, where <p a and <pp are the 
bulk water compositions in the two phases. This quan- 
tity is the so-called Gibbs transfer free energy per ion, 
which determines the ion partition and a Galvani poten- 
tial difference between the two phases. Its magnitude 
typically much exceeds fc^T (about 15fcsT for Na + and 
Cl _ in water- nitrobenzene at room temperatures). Fur- 
thermore, the dissociation of ionizable groups on the col- 
loid surface should be treated as a chemical reaction sen- 
sitively depending on the local environment (the solvent 
composition and the local elect ric p otential) as in the case 
of polyelectrolytes [llj, [l2j, |22j, |23j . Therefore, even very 
small composition variations around the colloid surface 
can induce significant changes in the ion distribution, the 
electric potential, and the degree of ionization, so it can 
drastically alter the colloid interaction. 

Historically, the colloid interaction has been supposed 
to consist of the screened Coulomb repulsive infraction 
-Pdlvo an d the van der Waals attractive interaction -F v dw 
since the celebrated theory developed by Derjaguin, Lan- 
dau, Verway, and Overbeek (DLVO) [lHj]. For two col- 
loid particles with radius a, the former reads 
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where d is the distance between two colloid centers, Q 
is the average charge on a colloid particle, e is the av- 
erage dielectric constant (at the average composition for 
a mixture solvent), and k is the Debye wave number. 
The latter arises from the pairwise van der Waals inter- 
action (oc — 1/r 6 ) among constituent molecules, where r 
is the distance between two molecules. In terms of the 
Hamaker constant An, -F v dw between two identical col- 
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loids with radius a is written as [3] 
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It grows for short separation £ — d — 2a <C a as 
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while it decays as — WA^a 6 /9d 6 for long separation I 3> 
a. We assume that F v d w (d) should saturate at a short 
distance on the order of the solvent molecular diameter do 
(~ 3A). Then F vdw (d) deer eases down to —Ana/12ao at 
contact. The size of An strongly depends on binary mix- 
tures and colloid particles under investigation and can be 
made small by matching the dielectric constants of the 
colloid and the solvent. In this paper, F v d w will not be 
included in our theoretical scheme for simplicity. 

In binary mixtures, the adsorption-induced composi- 
tion disturbances give rise to an attractive interaction 
F a( j between two colloid particles @. A linear theory 
can then be developed when the adsorption and the ion- 
ization are both weak. Further in the special case of 
weak selective solvation, the pair interaction is the sum 
of .Fdlvo in Eq.(1.2) and F a( j given by 



Fid — — ^adO £ 2 



2 d' 



(1.5) 



where £ is the correlation length growing near the sol- 
vent criticality. The coefficient A^d is proportional to 
h^, where h% is the surface field arising from preferen- 
tial molecular interactions between the surface and the 
two solvent species [25T4271 ] . In deriving Eq.(1.5), hi is 
assumed to be small and Eq.(1.5) is invalid very close to 
the criticality [Hj]; nevertheless, F a d becomes the dom- 
inant interaction at long distances for £ > We 
shall see that can exceed Frdw even at the molecular 
separation ao for small Ah- Furthermore, on approach- 
ing the solvent criticality without ions, the adsorption- 
induced interaction among solid objects (plates, rods, 
and spheres) becomes universal (independent of the ma- 
terial parameters) in the limit of strong adsorption 12 9| , 
so it has been called the critical Casimir interaction [3CH — 
l33l ]. However, it should be affected by ions for strong 
selective solvation. 

In this paper, we aim to investigate the ion effects on 
the colloid interaction in binary mixtures using a coarse- 
grained approach. A merit of our approach is that we 
can treat the preferential solvation in its strong cou- 
pling limit. In our recent work [l9l . |2~4| . we found that 
a strongly selective solute can induce formation of do- 
mains rich in the preferred component even far from the 
solvent coexistence curve. We shall see that this pre- 
cipitation phenomenon can occur on the colloid surface, 
leading to a wetting layer coating the colloid surface. As 
another prediction, there can be a first-order prewetting 
surface transition between weak and strong adsorption 
far from the solvent criticality, as discussed in our paper 



on charged rods 23j. These two phase transitions occur 
when the volume fraction of the selected component is 
relatively small (</> < <p c ). With intensified mutual in- 
teractions, colloid particles should trigger a macroscopic 
phase separation to form a floccuated phase [f| H, Q . As 
a closely related example, precipitation of DNA has been 
observed with addition of ethanol in water [33 - l36l | , where 
the ethanol added is excluded from condensed DNA. 

In addition to usual hydrophilic ions, we are also inter- 
ested in the colloid interaction in the presence of antag- 
onistic ion pairs in aqueous mixtures, where the cations 
are hydrophilic and the anions are hydrophobic, or vice 
versa [HI, UH . A well-known example in electrochemistry 
is a pair of Na~ and tetraphcnylborate BPhJ, where 
the latter anion consists of four phenyl rings bonded 
to an ionized boron and acquires strong hydrophobicity. 
Such ion pairs behave antagonistically in the presence 
of composition heterogeneities, giving rise to formation 
of mesophases, as recently observed by Sadakane et al. 
by adding a small amount of NaBPh.4 to D2O and tri- 
methylpyridine |37| . They should produce an oscillatory 
interaction between walls or colloid particles as in the 
case of liquid crystals [38| . 

The organization of this paper is as follows. In Sec. II, 
we will present a Ginzburg-Landau model of a binary 
mixture containing ions and ionizable colloid particles, 
where the bulk part includes the electrostatic and sol- 
vation interactions and the surface part the dissociation 
free energy. In Sec. Ill, we will examine the linearized ver- 
sion of our theory for the electrostatic and composition 
fluctuations as a generalization of the Debye-Huckel and 
DLVO theories. In Sec. IV, we will discuss how a wet- 
ting layer is formed on the colloid surface, which takes 
place as a precipitation phase transition. In Sec.V, we 
will present numerical results on the basis of our nonlin- 
ear scheme, where we shall encounter precipitation and 
prewetting phase transitions on the surface even far from 
the solvent coexistence curve. We shall also see bridging 
of wetting layers and a changeover between complete and 
partial wetting. 



II. GINZBURG-LANDAU MODEL FOR A 
MIXTURE SOLVENT 



We suppose monovalent hydrophilic cations and anions 
in a binary solvent composed of a water-like polar com- 
ponent (called water) and a less polar component (called 
oil) in a cell with a volume V [la, G3, H3| ■ We also place 
one or two negatively ionizable colloid particles in the 
cell. Experimentally, we suppose a dilute suspension of 
colloid particles in a mixture solvent. To apply our re- 
sults to such systems, we should set V — l/n co i, where 
n co i is the colloid density. 

For simplicity, we neglect the van der Waals interaction 
F v dw The Boltzmann constant ks will be set equal to 
unity in the following. 
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Free energy including electrostatics, solvation, 
and surface interaction 



We assume that the counterions coming from the col- 
loid surface are of the same species as the cations added 
as a salt. The cation and anion number densities are writ- 
ten as n\ and 712, respectively, while the water composi- 
tion is written as <j>. We treat these variables as smooth 
functions in space. The total free energy F to t — F + F s 
consists of the bulk part F and the surface part F s . The 
former is written as 

F = f dr[f tot + ^\Vcj ) \ 2 }+ J dr £ -^, (2.1) 

where J dr is the space integral in the colloid exterior in 
the cell and J dr is that in the whole cell including the 
colloid interior. The electric potential $ is defined even in 
the colloid interior. We assume that inhomogeneity of (j> 
gives rise to the gradient free energy, where the coefficient 
C is a positive constant. 

In the first term of Eq.(2.1), the free energy density 
/tot is the chemical part depending on 0, n\, and n% as 



/tot = /(<«+ T^, 



ln(niAf ) - 1 - gi4> 



(2.2) 



In this paper, the molecular volumes of the two solvent 
components take a common value vq, though they are of- 
ten very different in real binary mixtures. As a molecular 
length, we introduce 



(2.3) 



which is supposed to be of order 3 A. Then C ~ a^ 1 
[39l | . The colloid radius a is much larger than qq. (In our 
previous papers |16l4l9ll24| . a has been used to denote the 

1/3 

molecular length v ' .) We neglect the volume fractions 
of the ions assuming their small sizes. In our numerical 
analysis, we adopt the Bragg- Williams form |40jj, 



T 



f = In <t> + (1 - 0) ln(l - 0) + X 0(l - 0), (2.4) 



where x = xCO is the interaction parameter depending 
on T. The critical value of x is 2 without ions. In the 
ionic part of Eq.(2.2), A, = h(2Tr/m i T) 1 / 2 is the thermal 
de Broglie wavelength of the species i with rrii being its 
mass. The dimensionless parameters 01 and <?2 represent 
the degree of selective solvation [H|,[l7j], i n terms of which 
the solvation chemical potential is /i* ol = const. — Tgi<j> 
and the Gibbs transfer free energy is Tgi{<f> a — 4>p) for 
the ion species i (see Eq.(l.l)). If <fi is the water compo- 
sition, we have gi > for hydrophilic ions and gi < for 
hydrophobic ions. In many aqueous mixtures, the am- 
plitude |<7j| well exceeds 10 both for hydrophilic ions and 
hydrophobic solutes (ljjj . 

The last term in Eq.(2.1) is the electrostatic part [4llj . 
where E = — V$ is the electric field. The electric poten- 
tial $ is defined in the whole region including the colloid 



interior. We assume continuity of U through the colloid 
surface neglecting surface molecular polarization. It fol- 
lows the Poisson equation, 



V • D = -V ■ eV$ = 4ttp, 



(2.5) 



where D — eE is the electric flux density. The dielectric 
constant e(</>) is assumed to be a linear function of <p in 
the solvent. Thus it behaves as 



e{4>) = £0 +£1 



(colloid exterior) 
(colloid interior), 



(2.6) 



which is sq in oil (at = 0), £0 
and e p in the colloid interior, 
monovalent. Then the charge density p is written as 



- £1 in water (at = 1), 
Let all the charges be 



p — e(ni — U2) — ead s 



(2.7) 



The first bulk term is nonvanishing in the colloid exte- 
rior and the second part arises from the areal density a of 
the ionized groups on the colloid surface, where S s is the 
delta function nonvanishing only on the colloid surface. 
There is no charge density in the colloid interior. As a 
result, there arises a discontinuity in the normal compo- 
nent v • _D, where v is the outward normal unit vector 
on the colloid surface. Let D + and D_ be the values 
of D immediately outside and inside the colloid surface, 
respectively. Then, 



(D 



D 



-Airecr. 



(2.8) 



On the other hand, on the cell boundary, we assume no 
surface free energy and no surface charge. In our simu- 
lation, we thus set 



v h ■ V0 = 0, 

v h ■ E = -v h ■ V<I> 



0. 



(2.9) 



where v b is its normal vector of the cell boundary. 

The density of the ionizable groups on the colloid sur- 
face is written as <7o . The fraction of ionized groups or the 
degree of ionization a is defined in the range < a < 1. 
The density of the ionized groups is written as [231 ] 



a a- 



(2.10) 



We treat a as a fluctuating variable depending on the 
local composition and potential. The surface free energy 
F s depends on (f> and a as 

F s = J dS(Tj^ + f d ), (2.11) 

where J dS is the integration on the colloid surface. Here 
we neglect the second order contribution (oc </> 2 ) present 
in the original theory [25| to the surface free energy den- 
sity (though it is relevant near the critical point for neu- 
tral fluids [26]). The coefficient 7 represents the short- 
range interaction between the mixture solvent and the 
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colloid surface (per solvent molecule) [25j] . We call 7 the 
surface interaction parameter or the surface field (though 
h\ = — T7 is usually called the surface field in the liter- 
ature H^, H3)- The fd in Eq.(2.11) is the dissociation 
(ionization) free energy density of the form [IJ [2^, [H| , 



u_ 

Tar 



= a In a + (1 - a) ln(l - a) + a(A - Ai<£), (2.12) 



where the first two terms arise from the entropy of select- 
ing the ionized groups among the ionizable ones, while 
Ao — A±<f> is the composition-dependent ionization free 
energy divided by T. We suppose that the ionization is 
much enhanced with increasing the water content, which 
means that Ai should be considerably larger than unity. 



B. Equilibrium relations 

In our finite system, the cation number increases with 
an increase of ionization, while the numbers of the an- 
ions and the solvent particles are fixed. Let no be the 
average density of the added salt and </> be the average 
water composition. They are important parameters in 
our problem as well as x- Then, 



dr(m(r) — no) — / dSa, 



dr{n 2 {r) - n ) = 0, 



dr(cj)(r) - (f>) = 0. 



(2.13) 



(2.14) 



(2.15) 



The right hand side of Eq.(2.13) is the number of the 
counterions from the colloid surface. These relations are 
consistent with the expression for the charge density p in 
Eq.(2.7). In equilibrium we should minimize the grand 
potential f2 defined by 



F- J drihtjy + Y^ 



(2.16) 



Here h, and \xi are introduced as Lagrange mul- 
tipliers owing to the constraints (2.13)-(2.15). They 
have the meaning of the chemical potentials expressed 
as h — 5F/5<p and fa = SF/Srii. 

To minimize f2, we superimpose infinitesimal devia- 
tions 5(f), Sni, 5n,2, and 6a on <f>, m, ri2, and a, respec- 
tively. First, we calculate the infinitesimal variation of 
the electrostatic part F e = J dreE 2 /8tt in F in Eq.(2.1). 
From the relation S(sE 2 ) = 2E ■ SD - E 2 Se, we obtain 



SF R = 



dr 



E 

STT 



dS eaoda 



dS <Pv b ■ {eSE + 5s E), 



(2.17) 



where the integration is in the colloid exterior in the first 
term, on the colloid surface in the second term, and on 



the cell boundary in the third term. From Eq.(2.9) we 
have Vb ■ E = Ub ■ 5E = on the collicd surface, so the 
third term vanishes. From Eqs.(2.1) and (2.17) we obtain 

h= ficfy-TCV^-TY.gm-^E 2 , (2.18) 



Hi = T]n(n i \l)-Tg i( j> + Z i e$, 



(2.19) 



where /' = df/dcj) in Eq.(2.18) and Z x = 1 and Z 2 = -1 
in Eq.(2.19). It folows the modified Poisson-Boltzmann 
relations for the ion densities, 



(2.20) 



where n® — X i 3 exp(/ij/T) are constants. We introduce 
the normalized electrostatic potential, 



U = e$/T. 



(2.21) 



Using the above h and /ii , we calculate the incremental 
changes of F and F s as 



5F = J dr[h5(p + ^ ^Sn,} 

i 

dS[{Cu ■ V(f>)6<p + {e$a )6a], 



5F S =T J dS'(7-crAi)(50 



+Ta / dS 



In- 



1 



A - A x < 



5a. 



(2.22) 



(2.23) 



Vanishing of the surface terms proportional to 5<f> in 5F + 
5F S yields the boundary condition of cf> on the colloid 
surface written as [2J| 



Cv ■ Vcf) = 7 — Ai<7oa, 



(2.24) 



where v is the outward normal unit vector on the col- 
loid surface. In the same manner, vanishing of the terms 
proportional to 5a yields [23| 



1 - a 



exp(-A + A 10 + U- m/T). 



(2.25) 



We multiply Eq.(2.25) by n\ in Eq.(2.20) at the surface 
to obtain the mass action law on the surface, 



arii 



(2.26) 



where the factor exp(C/ — fi\/T) is cancelled. We intro- 
duce the composition-dependent ionization constant by 



K{cj>) = Ar 3 exp[-A + (A, 



9i) 



(2.27) 



in terms of which we have a = K(cf>)/[ni + K(cf>)]. Thus 
we have weak ionization a <C 1 for n\ 3> K ((f)) and strong 
ionization a = 1 for ri\ -C K((f>) on the surface. 

In addition, from Eq.(2.25), the ionization free energy 
density fd in Eq.(2.12) becomes 



fd/Tao = a{U - //i/T) + ln(l - a), 
which will be used in deriving Eq.(3.17). 



(2.28) 
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C. Changeover from hydrophobic to hydrophilic 
surface with progress of ionization 

We further discuss the consequence of the boundary 
condition (2.24). Near the surface, oil is enriched for 
<j>' = v ■ V0 > 0, while water is enriched for cf>' < 0. For 
very small a, the colloid surface is hydrophobic for 7 > 
and is hydrophilic for 7 < 0. However, with increasing a, 
an originally hydrophobic surface can become effectively 
hydrophilic if 



in terms of <j> and uq. From Eq.(2.20) we obtain 



< 7 < Aicr , 



(2.29) 



under which the surface derivative </>' becomes negative 
for a > 7/A1CT0. In contrast, if 7 > Aitro, the surface 
remains hydrophobic even for a = 1. 

This weakly hydrophobic situation can well happen in 
real colloid systems in mixture solvents for not small 
(To owing to strong composition-dependent ionization. 
As stated in Sec. I, colloid aggregation in near-critical 
lutidine-water occurred at lutidine-rich compositions for 
small ionization and at water-rich compositions for larger 
ionization 0, 0] ■ In our theory this means that the col- 
loid surface remained hydrophobic for small ionization, 
while it became hydrophilic for large ionization. It is 
also well-known that hydrophobic polyelectrolytes (with- 
out ionization) can become hydrophilic with progress of 
ionization jl3|. 



III. LINEAR THEORY IN ONE-PHASE STATES 

The Debye-Hiickel and DLVO theories [jHiJ are justi- 
fied for small electrostatic perturbations, where the am- 
plitude of the normalized potential U = e$/T should be 
smaller than unity. Here we present a generalized linear 
theory, including the composition fluctuations in a mix- 
ture solvent. To justify the linear treatment, we assume 
that the degree of ionization a and the surface field 7 are 
both very small. Treating a and 7 as small parameters, 
we calculate U and the composition deviation, 



ip = cf> — <f>. 



(3.1) 



to linear order in a or 7. From Eq.(2.24) the preferred 
component is only weakly adsorbed on the colloid surface. 
These deviations produce a change in the grand potential 
Q in Eq.(2.16) of second order (oc a 2 , a-f, or -f 2 ). 



A. Linearized relations 

In the limit of large cell volume V ^> a 3 or in the dilute 
limit of colloid suspension, we may assume that <f>, rii, 
and U tend to <fi, hq, and 0, respectively, exponentially 
far from the colloids. From Eqs.(2.18) and (2.19) we then 
have h = f'{(j>) -Tj^i ff<Mo and m = Tln(n Af ) - Tg^>. 



8n\ 
Sn 2 



m 

n 2 



no{giip 
n (g2ip ■ 



U), 
U). 



(3.2) 



From Eq.(2.25) the expansion of a is of the form, 

a = a[l + Aii/j + U], (3.3) 

where the surface values of tp and U are used and a — 
K(<fi)/no is the degree of ionization in the homogeneous 
case assumed to be very small. The deviation 5a = a — a 
is already of second order. We shall see that 8a gives rise 
to a third order contribution to Q and may be neglected 
in our linear theory. The colloid surface is under the 
fixed charge condition in the linear theory. We should 
calculate U in the whole space imposing V 2 U — in the 
colloid interior, while tp, n\, and n 2 are defined only in 
the colloid exterior. 

In terms of the average dielectric constant e — e((f>), 
we introduce the Bjerrum length £b and the Debye wave 
number k by 



i B = e 2 /eT, 

k = {8im Q e 2 /£T) 1/2 



(8^ B n ) 1/2 . 



(3.4) 
(3-5) 



Without coupling to the ion densities, the correlation 
length of the composition fluctuations is given by £ = 
(C/r) 1 / 2 in one-phase states with 



T 



1 



1 



wo 0(1 - <t>) 



2xL 



(3.6) 



where f"{4>) = d 2 f/d(j) 2 and use is made of Eq.(2.4). 

In the bulk region of the colloid exterior, Eqs.(2.5) and 
(2.18) are linearized with respect to tp and U as 

V 2 U = k 2 (U -g a ip), (3.7) 
CV 2 ^ = (r - r c )i, + 2n Q g a {U ~ g a ^), (3.8) 

where we introduce two coefficients, 

t c = n (gi +.g 2 ) 2 /2, (3.9) 

fja = (91 ~ 00/2. (3.10) 

As can be seen from the structure factor of the composi- 
tion in Appendix A, r c is the shift of the spinodal in the 
long wavelength limit due to the selective solvation. The 
size of t c can be significant even for small no for gx ^> 1 
and <?2 ^* 1) which agrees with experimental large shits 
of the coexistence curve induced by hydrophilic ions [42| . 
We assume r > r c to ensure the thermodynamic stability. 
The correlation length of ip is changed by ions as 



-1/2 



(3.11) 



which grows as r — > t c or as \ ~> Xs = 1/20(1 — (j>) — 
vqt c /2. The g a arises from the asymmetry of the selective 
solvation between the cations and the anions, giving rise 
to the coupling of U and ip. 
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On the colloid surface, Eqs.(2.8) and (2.24) yield the 
linearized boundary conditions, 

e{u ■ VU)+ - e v (y ■ VJ7)_ = Ane 2 a/T, (3.12) 
(v ■ V^)+ = 7/C, (3.13) 

where (•••)+ and (•■•)- denote taking the values im- 
mediately outside and inside the colloid surface, respec- 
tively. The a and 7 are averages defined by 



a = aoa, 

7 = 7 — A\<J. 



(3.14) 
(3.15) 



In fd in Eq.(2.28) we use the expansion ln(l — a) = —a + 
a 2 /2 + ■ ■ ■ , where a behaves as in Eq.(3.3) and a 2 = a 2 . 
Up to the second order we find 

fd/Tao = -a/ii/T -a- aAnp + a 2 /2. (3.16) 

Thus the last two terms in the grand potential in 
Eq.(2.16) become F s + J dS^\a = T J dSjip+const. in 
terms of 7, where the first term in the right hand side 
of Eq.(3.16) cancels to vanish. Hence the second-order 
contributions to arc written as 



+ JdreE 2 /8nT + JdSjip 



dSfrip - aU), 



(3.17) 



where the bulk integrations in the first two lines are trans- 
formed into the surface ones in the third line with the aid 
of Eqs.(3.7) and (3.8). With the third line, we thus need 
to calculate only the surface averages of t/j and U at fixed 
surface charge in the linear theory. 



B. Two characteristic wave numbers q\ and q2 

In the colloid exterior, there arise two characteristic 
wave numbers, denoted by qi and q 2 - If 9a = 0, U varies 
on the scale of the Debye screening length where k 
is defined in Eq.(3.5), and tp varies on the scale of £ in 
Eq.(3.11). For g a ^ 0, they are expressed as 



,1/2 ,1/2 



(3.18) 



From Eqs.(3.7) and (3.8) Ai and Ai are the solutions of 
the quadratic equation, 

A 2 - (At 2 + 1 - 7p)A + M 2 = 0. (3.19) 

Therefore, Ai and Ai satisfy 

Ai + A 2 = M 2 + 1 — 7p, AiA 2 = M 2 . (3.20) 



We also find (Ai — 1)(A 2 — 1) = J 2 , which will be used in 
the following calculations. Here the two parameters M 
and 7 P are defined by 



M = lM=[(r-r c )/C] 1 / 2 /^ 
7 P = | 3a |/(47r^ B ) 1/2 , 



(3.21) 
(3.22) 



where M — > as r — > t c and 7 P conveniently represents 
the solvation asymmetry of the two ion species. That 
is, 7 P should be smaller than unity for usual hydrophilic 
ion pairs [H, [ltl In Appendix A, we shall see that 
our system is linearly unstable for 7 P > M + 1 against 
charge-density-wave formation, so we limit ourselves to 
the region 7 P < M + 1 in one-phase states. In Fig.l, we 
display q\/n and qij^ in the M-"f p plane. 

Some typical cases are as follows, (i) As g a — > 0, we 
have q% — K and (72 — £ < (ii) Near the criticality, 
we have > 1 or M < 1. Furthermore, supposing 
hydrophilic ions, we assume that 7 P is not close to unity 
and the inequality 1 — 7 p ^> M holds. We then find 



91 - K (! -7 P ) 



2N1/2 



<12 



r x (i 



-1/2 



(3.23) 



where q\ ^> (72 or Ai 3> A2. (in) As shown in Fig.l, q\ 
and q2 are complex conjugates in the region \M — 1 < 
7 P < M + 1. From Eq.(3.19) the real part q R =Re(gi) 
and the imaginary part qi =Im(qi) are calculated as 



(3.24) 



te = [(A/ + l) 2 - 7p 2 ] 1 /2 K /2 ; 

9/ = [7p 2 -(M-1) 2 ] 1 / 2 ^/2. 

Oscillatory behavior appears for qi > q R or for 7 P 2 > 
M 2 + 1. As we approach the spinodal line or as M — > 
7 P — 1 with 7 P > 1, qn becomes small as (M — 7 P + l) 1 ' 2 , 
while qi tends to nonvanishing (n/t^) 1 / 2 . 

C. Profiles around a single colloid particle 

In the presence of a single colloid with radius a, we ob- 
tain the fundamental profiles U — Uo(r) and ip = ipoir) 
induced by nonvanishing 7 and a from the boundary 
conditions (3.12) and (3.13). In the large system limit 
F » a 3 , they are expressed as linear combinations of two 
Yukawa functions, e~ qi ^ r ~°^ jr and e~ q2 ^ r ~°^ /r, where r 
is the distance from the colloid center. As will be shown 
in Appendix B, they are of the forms, 



Un = 



g a Bie 



-q x (r-a) 



9aB 2 e 



-q 2 (r-a) 



(1 -Ai)(l + gi a)r (1 - A 2 )(l + q 2 a)r ' 

B ie -9i(r-o) S 2e -92(r-a) 



(1 + qia)r (1 + q2a)r 
The coefficients Bx and B 2 are defined by 
B t = a 2 {j(l - \ t ) - ag^/C^ - \ 2 ). 



(3.25) 



(3.26) 



In Appendix A, we will also express the correlation func- 
tions of the composition and the ion densities as linear 
combinations of these Yukawa functions. 
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FIG. 1: (Color online) Real and imaginary parts of qi/n and 
qz/n in the plane of 7 P and M = in the linear theory. 

Here qi and q2 are positive for 7 P < \M — 1] (in blue) and 
are complex conjugates for \M — 1| < 7 P < M + 1 (in red), 
while the system is unstable for 7 P > M + 1 (in white). These 
three regions are depicted in the right bottom panel, where 
the profiles and the interaction free energy will be shown at 
four x points (A), (B), (C) and (D) in Fig.2 and those at 
three • points (E), (F), and (G) in Fig. 3 in the linear theory. 



We examine some limiting cases, (i) As g a — ► 0, we 
find Ai — > 1, A 2 — > M 2 , and 



U = Q 



eBe - K (r-a) 

(1 + na)r 



ip = -7 



a 2 e -(r-a)/i 

c(i+ a /iy 



(3.27) 



where Uq is the Debye-Hiickel form [H-Q with Q being 
the average charge of a colloid particle, 



Q 



-Ana ae. 



(3.28) 



(ii) We assume \qt\a <C 1 and \qi\(r — a) <C 1, which hold 
in the limit of small Even for g a ^ 0, Eq.(3.25) yields 
the Coulombic behavior, 



U = Qt B -, Vo = -(ja 2 /C)-, (3.29) 



which follow from the relations B\ — B 2 = —7a 2 /C and 
Bi/(1 - Ai) - B 2 /(l - A 2 ) = -Anl B oa 2 jg a . (iii) Let us 
approach the instability line M = 7 P — 1 in Fig.l. In 
this case, both Uq and ip grow as Im[(Ai — A 2 ) _1 ] = 
— \/2qiqji oc (M — 7 P + l) -1 / 2 . In this limit, the linear 
theory is valid only for very small a and 7. 



D. Interaction between two colloid particles 

We suppose two colloid particles of the same species 
with radius a. In Appendix B, we will derive the inter- 
action free energy i*] n t from Afi in Eq.(3.17) as a linear 
combination of two Yukawa functions, e~ Q1 ( d ~ 2a ' / d and 
e -qi{d-ia) 1^ wnere d j s the separation distance between 
the two colloid centers longer than 2a. We express it as 



-Pint 



47ra 4 



(l + qia) 2 d (l + q 2 a) 2 d 
where the coefficients E\ and E 2 are defined by 



[7(1 - Aj) - ag a } 2 
C(l-Ai)(Ai-A 2 )' 



(3.30) 



(3.31) 



Notice that Fj nt is independent of the colloid dielectric 
constant e p . In Appendix B, we can see that e p appears 
in the third order contribution. 

Some limiting cases are as follows, (i) In the limit of 
weak solvation g\ — > and g 2 — > 0, we have E\ — >■ A-kIbct 2 
and E 2 — > 7 2 /C, leading to a decoupled expression, 



Fiat = -Fdlvo + Fad 

Q2 e -K(d~2a) ^ Ta i- 2e -(d~2a)/d 



6(1 



47T- 



C(l + a/£) 2 d 



(3.32) 



The first term is the DLVO interaction i^DLVO in Eq.(1.2) 
- The second term represents the adsorption- 
induced attraction F a d for small 7. For neutral colloids, 
we have 7 = 7 and h\ = — T7 to find the expression (1.5) 
with j4 a d = 4irTj 2 /C. Note that the linear theory is not 
applicable very close to the criticality [28[, as stated be- 
low Eq.(1.5). (ii) When \qi\a <C 1 and \qi\{d - 2a) <C 1, 
we use E\ — E 2 = Ai:Ibo 2 — 7 2 /C to obtain 



F int = AnTa\\nl B a 2 - f/C)/d, 



(3.33) 



(iii) Near the criticality and for hydrophilic ions, we may 
assume M <C 1 and 1 — 7 p 3> M, where qi and q 2 are 
given by Eq.(3.23). Then Fj nt in Eq.(3.30) takes the same 
form as the decoupled expression (3.32): 



-Pint — 



Q2 e - qi (d-2a) 

e(l + qia) 2 d 



47T- 



Ta 4 j 2 e- q ^ d - 2a ^> 
C(l + q 2 a) 2 d ' 



(3.34) 



where Q and 7 in Eq.(3.32) have been replaced by Q e 
and 7 e defined by 



= Q(l-g a i/^C£ B a)/ x /l~ 1 2 , 



7e = (7 - °9a)l\j\ -7p- 



(3.35) 



(v) When gi and g 2 are complex conjugates in the region 
|M - 1| < 7 P < M + 1 in Fig.l, Eq.(3.30) gives 

F\at = 87rTa 4 exp[-q R (d - 2a)] - 

x{ Jcos[g/(d - 2a)] + i4Tsin[a 7 (d - 2a)]}, (3.36) 
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FIG. 2: (Color online) Results of the linear theory with 
a — 600ao for hydrophilic ions with j p — 0.23, g a = 2, 
a = O.OOlaQ 2 , and k = O.laQ 1 , where M is (A) 0.4, (B) 0.5, 
(C) 1, and (D) 3 (see Fig.l for their locations in the M-7 P 
plane). Top: Normalized potential Uo(r) (left) and compo- 
sition deviation ipo(r) (right) in Eq.(3.25) vs r/ao around a 
colloid. Left bottom: Interaction free energy Fi n t(d) between 
two colloids in Eq.(3.30) vs d/ao- Right bottom: Minimum of 
Fi n t, denoted by F££ n , in the M-'y/a plane with the common 
<7, k, and a. It is zero (attained at infinity) in the right (for 
M > 1.3) and is negative in the left (for smaller M). 



where the coefficients J and K are the real and imaginary 
parts of Ei/(l + qia) 2 . On approaching the spinodal line 
7 P = M + 1, K grows as 1/qn but J remains finite. 

For neutral colloids, we compare the van der Waals in- 
teraction F v d w (d) in Eq.(1.2) and the adsorption-induced 
interaction F a d(d) in Eq.(1.5) or (3.32) at the closest sep- 
aration d — 2a — do, where 7 = 7 and £ = £. If C ~ Oq 1 
13911 . we estimate their ratio as 



^vdw 



W 2 (ao ai ) 2 T/AH {a » £) 
W 2 (a fr) 2 T/A H (a«0. 



(3.37) 



Even at the closest separation, the van der Waals attrac- 
tion is negligible when Ah/T is smaller than (10aoa7) 2 
for a 3> £ and than (10ao£7) 2 f° r a ^ £• 



E. Plotting analytic results in the linear theory 

On the basis of the analytic expressions (3.25) and 
(3.30), we plot Uo(r) and V , o( r ) around a single colloid 




700 

r/a 



1200 1250 1300 1350 1400 1450 

d/a 



FIG. 3: (Color online) Results of the linear theory with 
a — 600ao for antagonistic ions with 7 P = 1.727, g a = 15, 
j/a = -1, and k = O.laQ 1 , where M is (E) 0.769, (F) 0.833, 
and (G) 0.909 (see Fig.l for their locations). Left: Normal- 
ized potential Uo(r)/a 2 a vs r/ao around a colloid. Right: 
Normalized interaction free energy Fi n t(d)/a 4 a 2 between two 
colloids from Eq.(3.30) vs d/ao. 



and -Fint between two colloids for various M. Assuming a 
large radius a = 600ao, we set £b = 3ao and k = O.la^" 1 , 

where ao = v^ 3 . Then no = 1.33 x lO -4 ^ 1 . 

In Fig. 2, with hydrophilic ions, we set 7 P = 0.23, 
g a = 2, a = O.OOlag 2 , and j/a = —2. We can see that 
ipo increases considerably with decreasing M, while Uq is 
rather insensitive to M . Remarkably, the curves of F{ nt vs 
d/a exhibit a negative minimum at an intermediate d for 
small M. With these selected parameters, the DLVO in- 
teraction Fdlvo is equal to 1.061 x lO^Ta 4 ^ 2 = 13.75T 
at d = 2ao- 

In Fig. 3, with antagonistic ions 0, [H, we dis- 
play U (r)/a 2 a and F int (d)/Ta 4 a 2 by setting 7p = 1.727, 
g a = 15, and j/a = -1. As M -> M c = 0.727, U (r) 
and F lnt (d) grow as stated below Eqs.(3.24) and (3.36). 
Oscillatory relaxation is conspicuous close to the instabil- 
ity. Here, the linear conditions, \Uo\ <SC 1 and \ipo\ <C 1, 
are ensured only by very small a and 7. Otherwise, the 
nonlinear theory is required. In Fig. 3, we thus divide Uq 
and i^nt by a 2 a and a 4 a 2 , respectively. 



IV. PRECIPITATION DUE TO HIGHLY 
SELECTIVE SOLVATION 

In a binary mixture in one-phase states, a highly se- 
lective solute can induce precipitation of domains rich 
in the preferred component [19|, |2J] • The solute can be 
either a hydrophilic salt (such as NaCl) or a neutral hy- 
drophobic solute. The equilibrium size of a precipitated 
domain depends on the system size V (as can be known 
from Eq.(4.8) below). 

In the next section, we will show numerically that pre- 
cipitation can occur on the colloid surface. Therefore, in 
this section, we summarize and extend our previous the- 
ory [24| to understand precipitation on the colloid sur- 
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face. We suppose an experiment of very dilute colloid 
suspension; then, the volume V assigned to each colloid 
particle is the inverse droplet density . For high col- 
loid concentrations, colloids should interact collectively 
due to the growth of wetting layers to form a floccuated 
phase I,!,©] . 

A. Bulk precipitation 

For hydrophilic cations and anions, the phase behavior 
of precipitation is little affected by charge density varia- 
tions, because they are significant only at interfaces, col- 
loid surfaces, or a container (see the right bottom panel of 
Fig. 10 as an example). Thus, neglecting the electrostatic 
interaction, we may use results of a three component sys- 
tem [24j by setting 

n 1 =n 2 = n/2, (4.1) 
9 = {91+92)11. (4.2) 

We assume the strong solvation condition j>1. In the 
precipitated phase, called a, the water composition <f> a is 
close to unity and the solute density n a is much larger 
than the average n = 2no by the factor e 9 '- 1- ^-' The 
precipitation effect is significant for not small g(l — </>). 

Let us decrease n in the presence of precipitated do- 
mains. Then the volume fraction 7 Q of phase a decreases 
and eventually tends to zero as 

7 „ = («K-1)^ (1 ^, (4.3) 

where n p is a minimum solute density for precipitation. 
Its asymptotic expression for g 3> 1 is given by 

n p = e-^-^G^,x)/T. (4.4) 
The function G(4>, x) is a positive quantity defined by 

Gfax) = -/(£> - (i - $)f\4>) 

= -(T/v )[\n^ + X (l-$)% (4.5) 

where the second line follows from Eq.(2.4). Alterna- 
tively, we may decrease x at fixed <j> and n in the presence 
of precipitated domains. The volume fraction 7 Q tends 
to zero as x approaches a lower bound Xp — Xpi^i^)- 
From Eq.(4.4), X p satisfies n = e~ 3 (1 ~^G(0, x p )/T. Use 
of the second line of Eq.(4.5) gives 

Xp = -[l n + e 9 ^-^v h]/(l - 0) 2 . (4.6) 

It also follows the relation, 

x - Xp = v Q e^\n - n p )/(l - 0) 2 . (4.7) 

Here n appears in the combination e 9 ^ 1- ^^. 

In Fig. 4, we display Xp vs 4> at no — 3 x lO -4 ^ 1 . 
The left panel gives the curves for g = 7,8,8.5,9, 9.5, 
and 10 in the range < < 0.47. With increasing g, 




0.1 0.2_ 0.3 0.4 ' 0.4 0.45 0.5 0.55 0.6 0.65 0.7 0.75 



<t> 4> 

FIG. 4: (Color online) Minimum interaction parameter \p vs 
<j> at n = 3 x lO^ 4 ^ 1 in the range < <p < 0.47 for g = 
7, 8, 8.5, 9, 9.5, and 10 from above (left) and in the range 0.4 < 
4> < 0.75 for g = 10 (right). Two-phase coexistence region (in 
yellow) without solute is in the upper part. Spinodal curve 
r = t c is also written by dotted line (right). 



the precipitation branch is suddenly detached downward 
from the solvent coexistence curve. For 5 > 1, they 
are in excellent agreement with the asymptotic formula 
(4.6) for Xp- (The asymptotic formula (4 4) for n p is 
also a good approximation for g 3> 1 Qji UM-) The 
right panel gives the curve for g = 10 at the same no 
in the range 0.4 < </> < 0.75. The curve is only slightly 
outside the coexistence curve for (f> > 0.6 and touches the 
spinodal curve r = r c at ( X ,$) = (0.5187,1.9728) 0. 
We obtained these curves numerically from minimization 
of the bulk free energy at (4>) = 4> and (n) — n, neglecting 
the surface free energy. Thus these curves are those for 
two-phase coexistence with a planar interface separating 
the two phases. 



B. Surface tension effect on a precipitated droplet 

The surface tension a s between a precipitated droplet 
and the surrounding solute-poor region is well-defined, 
though the precipitation is a nonlinear effect of a highly 
selective solute. It is of order a s ~ T/a^ far from the 
solvent criticality, being independent of the solute den- 
sity. We here examine the surface tension effect on the 
droplet stability. 

We consider a single spherical droplet of phase a with 
radius R in a large volume V. The droplet volume frac- 
tion is then j a = v/V, where v = AttR 3 /3. For n > n p 
and at very small volume fraction 7™ <C e _s< - 1_ ^\ the 
droplet free energy is expressed as 0, H3] 

\ pi 2V 
= Ana s (R 2 - 2R 3 /3R C + R 6 /3i&). (4.8) 

In the first line, the first two terms constitute the stan- 
dard droplet free energy in the nucleation theory and the 
third term (oc R 6 /V) arises from the finite size effect [40| . 
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In the second line, we introduce a critical radius R c and 
a minimum radius R m by 



R c = 2a s e- g{1 -^/T(n- 
= 2v Q a s /[T(l-4>) 2 (x 
R m = (3cr s V/2TTTn) 



Xp)}, 



l/4 e -9(l-4>)/ 2 . 



(4.9) 
(4.10) 



where Eq.(4.7) is used in the second line of Eq.(4.9), R c 
grows as x - > Xp> an d -Rm decreases with decreasing 4>. 
We estimate R m ~ (V"ao/-uo™) 1/4 e~ 9(1 ~* )/2 far from the 
solvent criticality. Minimization of the second line of 
Eq.(4.8) yields 



(R/Rm) 4 = R/R c - 1 > 



(4.11) 



for the equilibrium radius R. We require AF '/ ' A-k<j s R 2 = 
(2 - R/R c )/Z < to find R > 2R C . Hence R > R m from 
Eq.(4.11) and R m is the minimum radius of equilibrium 
droplets in a finite system [46j]. For R m 3> R c , droplets 
with radii much larger than i? m can appear as 



R 3 = R m / R c 



(4.12) 



where the right hand side is proportional to V and is 
independent of a s in agreement with Eq.(4.3). 

It also follows the condition R m > 2R C from Eq.(4.11) 



48] , leading to lower bounds of n 



n p > Aa s e~ 9 ^/TR m , (4.13) 
X-Xp >Av o a s /TR m (l-0) 2 , (4.14) 

for the formation of a droplet. Thus the bulk precipi- 
tation curves n — n p and X = Xp are shifted upward 
by amounts proportional to a s /R m oc o 3 J 4 V~ x l^ for 
droplets with surface tension. 

As an example, let us set V — Ax 10 9 7raQ/3 with C = 
2aQ 1 , g = 10, n = 6 x lO" 4 ^ 1 , and 4> = 0.41. The 
surface tension is then a s = 0.06T/O 2 , as \ ~ ► = 1-871 

M 



and x - Xp as 



Thus we obtain n p = 6 x 10 4 and R m — 35.1ao- 
The right hand side of Eq.(4.13) is 1.9 x 10~ 5 , while that 



of Eq.(4.14) is 0.02. 



surface and the wetting layer as a constant. As a general- 
ization of Eq.(4.8), R is determined by minimization of a 
free energy contribution _F we t- In terms of R c in Eq.(4.9) 
and R m in Eq.(4.10), it is expressed as 



iv 



Aira s 



R 2 



R„ 



2 3 
3R C 



1 



3Ri 



(R 3 



.3\2 



(q + s 3 ) 2 / 3 



2Rn 
3R r 



<t_ 
3 



(4.15) 



The first line tends to Eq.(4. 
line, q and s are defined by 



as a — > 0. In the second 



(R 3 - a 3 )/R. 
a/R m . 



3 

m; 



(4.16) 
(4.17) 



We treat q as an order parameter. In Eq.(4.15), the 
selective solvation is accounted for in R c and R m , but 
the electrostatic interaction is neglected. See Fig. 18 in 
Ref.[IJ] for q in the s-R m /R c plane (where s in Eq.(4.17) 
is written as D). In the thin layer limit R 3 /a 3 - 1 < 1 or 
for q <^ s 3 , F wct is expanded up to the third order with 
respect to q as 



4na s Rl V R c ' 3s K 



1 ^9 



Aq 3 



)— + 
' 9 81s 7 



(4.18) 



Here the coefficients of the first two terms can vanish 
at a/R c = 1 and s = 3 -1 / 4 , where we predict tricriti- 
cal behavior with varying x or Rc- That is, q changes 
continuously or discontinuously depending on whether 
a > i- 1/i R m or a < 3-^ 4 R m . 

For a > 3 _1 / 4 i? m , q becomes nonvanishing for 1 — 
a/R c < or for i?" 1 > a -1 continuously as a second- 
order phase transition. From Eq.(4.9) this condition 
yields a lower bound of x ~ Xp m the form, 



X ^ X P > 2v (T s /Ta(l 



(4.19) 



which also follows if R m is replaced by 2a in Eq.(4.14). 
With further increasing a much above i? m or for s 1, 
the first two terms in Eq.(4.18) gives q = (a/R c — l)/s or 



C. Wetting layer formation on a colloid surface 

A completely wetting layer can appear on a colloid sur- 
face above a precipitation curve for the hydrophilic case 
7 < or for the hydrophobic case 7 > under the condi- 
tion (2.29). In our numerical analysis, the precipitation 
curve for a colloid is only slightly shifted upward from the 
bulk curve x = Xp in the x~4> plane. We also recognize 
that the precipitation curve is nearly independent of the 
surface parameters 7 and 00 (see Figs. 6 and 8). 

We suppose that a colloid with radius a is com- 
pletely wetted by a spherically symmetric layer with 
thickness R — a. The volume fraction of phase a is 
7 Q = 47r(_R 3 - a 3 )/3V in a volume V. For R - a > a , 
we may treat the surface free energy between the colloid 



R 3 /a 3 - 1 = (R m /a) 4 (a/R c - 1), 



(4.20) 



which holds for a/R c — 1 (a/7? m ) 4 because we have 
assumed q <C s 3 in setting up Eq.(4.18). The situation 
a 3> Rm can well happen in real systems particularly 
for relatively small <\> (where R m is decreased). In such 
cases, the factor (i? m /a) 4 in Eq.(4.20) is small and the 
wetting layer thickens slowly with increasing x- On the 
other hand, in the thick layer limit R 3> a, R is given by 
Eq.(4.12). 

For a < 3 _1 ^ 4 i? m , q becomes nonvanishing discontinu- 
ously as a first-order phase transition when 1 — a/R c is 
smaller than a small positive constant (oc (s — 3 -1 / 4 ) 2 
near the tricriticality) . As a result, the lower bound of 
X — Xp f° r precipitation is slightly smaller than the right 
hand side of Eq.(4.19). 
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In the above theory, we have examined the transition 
between weakly adsorbed states and completely wetted 
states. In the next section, however, we shall see that 
a hydrophobic surface can be partially wetted by a non- 
spherical water-rich layer at relatively small cj> and for 
not very large 7 (see Figs. 16 and 17). 



V. NUMERICAL RESULTS 

In this section, we present numerical results on the 
basis of our nonlinear theory in Sec. II with a single colloid 
or two colloids placed at the center of a large cell. The 
correlation length £ in Eq.(3.11) increases up to 14.1ao at 
X = 1.965 in Fig. 12, but it is only a few times longer than 
ao in the other examples far from the solvent criticality. 



Parameter values selected 



We set C = 2<2q , A = 18, A x = 12, i B = 
e 2 /eoksT — 3a , and E\ = Eq. We use large Ai, so the 
degree of ionization a sensitively depends on the surface 
composition from Eq.(2.25). Several values are assigned 
to the density of ionizable groups, surface with 7 > 0. 
Because of the numerical convenience, the colloid radius 
a is assumed to be rather small as 



Hydrophobic surface 



a = 15ao- 



(5.1) 



We also performed simulations with a = 25ao to obtain 
essentially the same results, though the corresponding 
figures are not shown. 

In Subsecs.VB-D, we treat hydrophilic ions with 171 = 
11 and g 2 — 9 at v n = 3 x 10~ 4 in Figs. 5-17. Re- 
call that we have introduced the Debye wavenumber k in 
Eq.(3.5) and the asymmetry parameter 7 P in Eq.(3.22) 
as functions of the average composition <f>. Here, we have 
K = 0.127ag 1 and 7 P = 0.135 for 4> = 0.4. In Subsec.VE, 
we treat antagonistic ion pairs with gx = —g 2 = 13 
at vqUq — 10~ 3 in Figs. 18 and 19, where we have 
k = O^a^ 1 and 7p = 1.76 for 4> = 0.4. 

In the case of a single colloid, our cell is a large sphere 
with radius 10 3 arj in the spherically symmetric geome- 
try. In the case of two colloids, it is a cylinder with 
radius 256ao and height 512arj in the axisymmetric ge- 
ometry. Since a = 15arj, the colloid volume fraction is 
3.375 x 10~ 6 for a single colloid and is 1.341 x 10~ 4 for 
two colloids. For hydrophilic ions and at <f> = 0.41, R m 
in Eq.(4.10) is 35.1arj in the single colloid case [43| and 
13.0arj in the two colloid case. In agreement with the dis- 
cussion in Subsec.IVC, the layer formation due to precip- 
itation is discontinuous for the single colloid case (where 
a < 3 _1 / 4 i? m ) but is continuous for the two colloid case 
(where a > ■$- 1 l i R m ). 
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FIG. 5: (Color online) if>(r) = <j>(r) — (f> and U(r) vs r/ao 
around a colloid from the linear and nonlinear theories below 
the precipitation line at vono = 3x 10~ 4 . Upper plates: Those 
for a hydrophobic surface with 7 = O.OSa^ 2 at ao = 0.0l8ag 2 , 
where at the surface is smaller than 1/gi and \U\ is at 
most of order unity so that the linear results fairly agree with 
the nonlinear results. Lower plates: Those for a hydrophilic 
surface with 7 = 0.2a^ 2 at small ao = 0.003ar7 2 . From ip ~ 
0.3 as r — > a, the cations are considerably accumulated near 
the surface and U becomes positive in the nonlinear theory, 
while U remains negative in the linear theory. 



B. Comparison of the linear and nonlinear theories 

We compare the profiles of ip(r) = <fi(r) — <f> and 
U(r) = e$(r)/T from the linear theory in Sec. Ill and 
those from the nonlinear theory in Sec. II below the bulk 
precipitation curve x < Xp- The linear theory is based 
on the assumptions (3.2) and (3.3) and is thus valid only 
f° r M ^ l/lffil an d \U\ < 1 near the surface. 

In the upper plates of Fig. 5, we suppose hydrophilic 
ions and a hydrophobic surface with 7 = 0.08a^" 2 . The 



other parameters are x = 1-908, 



0.425, and a 



0.018<2q . Here a — 0.25 from the nonlinear calculation, 
leading to a = 0.0045 and 7 = 0.026, which were then 
used in the linear calculation. In this case, gi<f> and U are 
both decreased near the surface and their amplitudes are 
not large compared to unity. Thus the linear results are 
in fair agreement with the nonlinear results. 

In the lower plates of Fig. 5, we suppose hydrophilic 
ions and a hydrophilic surface with 7 = — 0.2a ' 2 . The 

= 1.869, 4> = 0.412, and er = 
0.992 in the nonlinear calcula- 
-0.236a,7 2 . In 



other parameters are x 
0.003ao 2 . We have a = 



tion, which gives a — 0.00297et and 7 
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this case, <j> is a few times larger than l/p, ~ 0.1 near the 
surface, leading to considerable ion accumulation, which 
cannot be accounted for in the linear theory. Further- 
more, the cations are more enriched than the anions near 
the surface, leading to positive U in the nonlinear theory, 
while U remains negative in the linear theory. 

Comparison of the linear and nonlinear theories will 
also be made in the presence of two colloids at </> = 0.5. 
See explanations of Fig. 12 below. 

C. Prewetting and precipitation on a colloid 



Hydrophilic surface and hydrophilic ions 
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With hydrophilic ions, there appear two transition 
lines of prewetting and precipitation for each colloid in 
the 4>-x plane. They are located far below the solvent 
coexistence curve for strong selective solvation. The 
prewetting line sensitively depends on 7 and rro (as will be 
seen in Figs. 6 and 8 below). It starts from a point on the 
precipitation line ending at a surface critical point, across 
which there are discontinuities in the physical quanti- 
ties. The transition across the precipitation line for a 
colloid is first-order for the present parameters. We fur- 
ther confirmed that the precipitation line exhibits no ap- 
preciable dependence on 7 and <tq and approaches the 
bulk one x — Xp with increasing a for a < 3 _1 / 4 i? m . If 
a > 3 _1 / 4 i? m , the precipitation transition becomes con- 
tinuous and it is difficult to precisely determine the loca- 
tion of the precipitation line. 



1. Hydrophilic surface with 7 = — 0.2a ' 



In our previous work 24J , we examined one example of 



a hydrophilic colloid with 7 = — 0.2aQ~ in the presence of 
hydrophilic ions, where precipitation and prewetting are 
first induced on the colloid surface before in the bulk re- 
gion. Again with 7 = — 0.2oq 2 , the left panel of Fig. 6 dis- 
plays three examples of the prewetting line corresponding 
to a 2 <7o = 0, 0.003, and 0.005, around which a is nearly 
equal to unity and a = <7r> A first-order precipitation 
line is also shown (broken line), which is independent of 
<7o- The right panel of Fig. 6 shows the preferential ad- 
sorption r around a colloid as a function of x f° r several 
4> at a 2 ao = 003. Since <j>(r) tends to a constant 4>{L) far 
from the colloid, T is calculated from 



r = J dr[tf,(r) - 4>(L)], 



(5.2) 



where the integration is in the colloid exterior. With 
precipitation, we have T ~ 4ir(l — 4>){R 3 — a 3 )/3 for a 
single colloid. 

In Fig. 7, profiles of <p(r) and the total ion density 
nx{r)+n-z{r) are displayed at four points (A)-(D) marked 



in Fig. 6, where 7 



-0.2an 



= 0.412, and a = 



0.003a^. We can see a weakly discontinuous prewetting 
transition between (A) and (B) and a strongly discontin- 
uous precipitation transition between (C) and (D). The 



FIG. 6: (Color online) Prewetting and precipitation transi- 
tions on a hydrophilic colloid surface with a — 15ao and 
7 = — 0.2o,q 2 . Left: Colloidal precipitation line (red broken 
line) and prewetting lines for a^ao = 0, 0.003, and 0.005, 
around which a = 1 here. Right: Normalized adsorption 
r/(47ra 3 /3) vs X at a 2 Q a = 0.003 for (a) cj> = 0.42, (b) 0.4176, 
(c) 0.415, (d) 0.412, and (e) 0.41. No phase transition appears 
on path (a), the prewetting critical point is on path (b), and 
the transition is discontinuous for paths (c), (d), and (e). 



Prewetting and Precipitation : a 2 ,-/ = -0.2 a^ao = 0.003 = 0.412 




FIG. 7: (Color online) <j>(r) (left) and v (ni(r) + n 2 (r)) (right) 
at four points (A)-(D) marked in Fig. 6 for a hydrophilic col- 
loid surface with 7 = — 0.2a^" 2 and cj> — 0.412. There is a 
weakly discontinuous prewetting transition between (A) and 
(B) and a strongly discontinuous precipitation transition be- 
tween (C) and (D). 



value of x at this precipitation transition is 1.8945, while 
it is 1.8926 for d = 25oo with the other parameters held 
at the same values (not shown here). These values of % 
are slightly larger than Xp — 1-8795 at (f> — 0.412. 



2. Hydrophobic surface with 7 = 0.08a 

The prewetting behavior is more exaggerated for a 
weakly hydrophobic surface under Eq.(2.29) than for a 
hydrophilic surface. In Subsec.IIC, we have discussed 
the changeover from a hydrophobic to hydrophilic sur- 
face with progress of ionization. Here, we examine the 
prewetting and precipitation transitions in the weakly 
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Hydrophobic surface and hydrophilic ions 



Prewetting : a 2 j =0.08 a 2 a = 0.018 



2.05 



two phase coexistence without salt 



precipitation with salt 




FIG. 8: (Color online) Phase diagram of prewetting and pre- 
cipitation transitions on a hydrophobic colloid with Oq7 = 
0.08. Three prewetting lines, corresponding to a^ao = 0.015, 
0.018 and 0.021, are more extended than in the hydrophilic 
case in Fig. 6. Precipitation on the colloid surface occurs above 
the broken red line. See Fig. 9 for a and V on vertical paths 
(a),. . . , and (g). 



hydrophobic case with 7 = 0.08<2q in Figs. 8-17 and also 
in the strongly hydrophobic case with 7 = 0.2a^~ 2 in the 
lower plates in Fig. 12. 

In Fig. 8, we display three examples of the prewetting 
line corresponding to a 2 <ro — 0.015, 0.018, and 0.021, 
which arc pushed downward with increasing <7rj. They 
even bend downward and much extend outside the bulk 
precipitation line. The precipitation line for a colloid 
(broken line) is inside the bulk precipitation region and 
is independent of <jq (as in Fig. 6). 

In Fig. 9, we show the degree of ionization a and the ex- 
cess adsorption T defined in Eq.(5.2) as functions of \ for 
various (j>. In the upper plates, the discontinuities across 
the prewetting line increase with increasing x, where the 
critical point is located at the smallest \ on the line. (In 
sharp contrast, they decrease with increasing x m the 
hydrophilic case in Fig. 6.) In the lower plates, a and T 
are shown along paths, (e), (f), and (g), where (e) and (f) 
pass the prewetting and precipitation lines but path (g) 
passes the precipitation line only. In particular, on path 
(e), a changes from about 0.2 to values slightly smaller 
than unity across the prewetting line or between (A) and 
(B), while r changes greatly across the precipitation line 
or between (C) and (D). The r/(47ra 3 /3) is -0.0137 at 
(A), 0.417 at (B), 1.30 at (C), and 19.1 at (D). 

The jump of F between (C) and (D) is very large and 
can be explained by minimization of -F W ct in Eq.(4.15) 

In Fig. 10, the profiles of 4>{r), U(r), ni(r), and n\(r) — 
n 2 (r) are shown at points (A), (B), (C), and (D) on path 
(c) in the right bottom panel of Fig. 9. At (A), which is 
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FIG. 9: (Color online) Degree of ionization a (left) and nor- 
malized adsorption r/(47ra 3 /3) (right) vs x on a hydrophobic 
colloid with ao7 = 0.08 and ao°o = 0.018 along paths (a) 
4> = 0.46, (b) 0.455, (c) 0.45, (d) 0.44, (e) 0.425, (f) 0.415, 
and (g) 0.41 (see their locations in Fig. 8). Top: Behavior 
around the prewetting line at Oqoq = 0.018 (middle line in 
Fig. 8) along paths (a),(b), (c), and (d). Bottom: Large jumps 
across the precipitation line along the paths (e),(f), and (g) 
(dotted lines). On path (e), there is a prewetting transition 
between (A) and (B) and a precipitation transition between 
(C) and (D). Path (f) also passes the two transitions. Path 
(g) passes only the precipitation line. 



below the prewetting line, the colloid surface is still hy- 
drophobic with a negative surface value of $ — d(j)/dr, 
U (r) is negative, and the ions are weakly accumulated 
near the surface. At (B), which is slightly above the 
prewetting line, the surface is hydrophilic, U has a small 
maximum, and the cations and the anions are both en- 
riched near the surface. At (D) there is a thick wetting 
layer enriched with ions. In the right bottom panel, the 
charge accumulation near the surface and the electric 
double layer (in the inset) are shown, which are small 
because of small difference g\ — gi = 2 in this case. In 
passing, the value of x a t the precipitation transition is 
1.9329, while it is 1.9320 for d = 25ao (not shown here). 
These values are only slightly larger than x p = 1.9246. 
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Prewetting and Precipitation : a„7 = 0.08 a„(To = 0.018 = 0.425 
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FIG. 10: (Color online) cj>(r) (left top), Uir) (right top), 
von\(r) on a semi- logarithmic scale (left bottom), and 
vo(ni(r) — fi2(r)) (right bottom) on a hydrophobic colloid at 
four points (A),(B), (C), and (D) on path (e) in Fig. 9, where 
a§7 = 0.08, ao<7o = 0.018, and (j> = 0.425. A prewetting tran- 
sition occurs between (A) and (B) , while a precipitation tran- 
sition between (C) and (D). At (B), (C), and (D), the ions are 
enriched in the layer. In the inset in the right bottom plate, 
a small electric double layer is shown in the interface region 
around r = 40an at (D). 



D. Two colloids and their interaction free energy 
for hydrophobic surface and hydrophilic ions 

Placing two hydrophobic colloids along the z axis, we 
now calculate the profiles of <fi and U and the interaction 
free energy Fi nt (d) as a function of the distance d between 
the two colloid centers in Figs. 11-17. The profiles depend 
on r = (x 2 + y 2 ) 1 / 2 and z. The degree of ionization 
a = a(9) depends on the angle 9 with respect to the z 
axis. Its angle average is written as 



{a) = / d6sm9 a(9). 
Jo 



(5.3) 



The precipitation transition is continuous for the present 
parameters, as stated at the beginning of this section. 



'sive to attractive interaction 



1. Crossover from 



First, for ef> = 0.425, we examine the behavior around 
the prewetting transition with Oq7 = 0.08 and o^gq = 
0.018. In Fig.ll, the left panel gives F lnt vs d for 



X = 1.88 and 1.92. Remarkably, it is small and posi- 
tive for x = 1-88 but is negative and is much amplified 
for x = 1.92. The right panel displays the corresponding 
<f)(r,z) at d = 44a n = 2.93a, where (r/(47ra 3 /3), (a)) = 
(-0.0198,0.305) for X = 1-88 and (1.30,0.941) for X = 
1.92 with a big difference in T. The screening length k' 1 
in Eq.(3.5) is 7.94ao, while the correlation length £ in 
Eq.(3.11) is 2.71a for X = 1-88 and 3.23a for x = 1-92. 
The surface remains hydrophobic ((/>' > 0) for x = 1-88 
but becomes hydrophilic ((f)' < 0) for x = 1-92. 

On the other hand, at the critical composition <fi = 0.5, 
there is no prewetting transition, though the precipita- 
tion occurs for x > Xp = 1-9675, as can be seen in Fig. 4. 
The instability point t = r c occurs at x = 1-97 since 
v t c = 0.06 (see Eqs.(3.9) and (3.11)). Thus, at = 0.5, 
the physical quantities change continuously with vary- 
ing x below Xp- In Fig. 12, we show Fi n t vs d in the 
left and the corresponding <p(r, z) at d = 44ao = 2.93a 
in the right for <j> = 0.5. The screening length k _1 is 
8.14ao, while the correlation length £ is 1.14a for x = 1.2 
and 14.1ao for x = 1-965. (i) The upper panels are 
for the weakly hydrophobic case with a§7 = 0.08 and 
a 2 a = 0.018, where {T / [Ana 3 /3), (a)) = (0.063,0.65) 
for x = 1-2 and (3.15,0.98) for x = 1-965. The sur- 
face is weakly hydrophilic at x — 1-2 and is strongly 
hydrophilic at x = 1-965. (ii) The lower plates are 
for the strongly hydrophobic case with ag7 = 0.2 and 
a§cr = 0.01, where (T / {Ana 3 /3), (a)) = (-0.093,0.64) 
for x = 1-2 and (-1.31, 0.32) for x = 1-965. The adsorp- 
tion of the oil component is enhanced on approaching the 
criticality. 

In the weakly hydrophobic case in the upper plates in 
Fig. 12, the linear theory fairly holds at x = 1-2, while 
it breaks down at x — 1-965. In fact, F- m t(d) /ksT at 
d = 2a is 3.71 (linear theory) and 2.38 (nonlinear theory) 
for x = 1-2, while it is -22.4 (linear theory) and -3.44 
(nonlinear theory) for x — 1-965. The surface value of 
]p = $-4, is about 0.02 for x = 1.2 and 0.25 for x = 1-965, 
while U ~ — 1 both for these cases. 



2. Bridged and disconnected wetting layers 

Next, we examine the wetting layer behavior after the 
precipitation transition in Figs. 13-15 for <f> = 0.425. In 
Fig. 13, the profiles of 4>{r,z) and U(r,z) are presented 
above the precipitation line with the same parameter val- 
ues as in Fig.ll. Here we have thick wetting layers on the 
two colloids and they are bridged for d = 80ao = 5.33a 
(top plates) and are disconnected for 102ao = 6.8a (bot- 
tom plates). In Fig. 14, the interaction free energy Fi nt (d) 
is shown as a function of d, where the three curves cor- 
respond to x = 1-94, 1.95, and 1.96. At relatively short 
separation d = 44ao, for instance, T increases with in- 
creasing x as F/(47ra 3 /3) = 3.87, 7.20, and 23.3, respec- 
tively. We recognize that the wetting layers are bridged 
for relatively small d but are disconnected for large d, 
exhibiting hysteresis. The i 7 l nt (d)/T assumes negative 
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FIG. 11: (Color online) Normalized interaction free energy 
Fint(d)/T vs normalized separation distance d/ao (left) and 
(f>(r, z) at d = 44ao = 2.93a (right) for <f> = 0.425, where af^-y = 
0.08, and agfo = 0.018. The system is below a prewetting line 
at x — 1-88 and above it at \ — 1.92, where Fi n t changes its 
sign from positive to negative. The surface is hydrophobic 
with (\> = v ■ \7<p > (see Eq.(2.24)) at % = !- 88 : but is 
hydrophilic with < due to an increase in a at x — 1.92. 
The color represents <f> according to the color bar attached. 



values of order —100 while bridged, but it becomes very 
small once detached. For larger colloid radii, the value 
of Fi ni (d)/T should be much more amplified (oc a 3 ). In 
Fig. 15, the hysteretic transition behavior is illustrated 
between the bridged and disconnected sates. That is, in 
an interval of d, we find two linearly stable profiles, where 
one is metastable with a higher i*i n t(d). 

It is worth noting that Hopkins et al. [l4[ also found a 
bridging of adsorption layers of two approaching neutral 
colloid particles in a mixture solvent close to the coexis- 
tence curve. 
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FIG. 12: (Color online) Normalized interaction free energy 
jFjnt(d) /UbT vs normalized separation distance d/ao (left) and 
4>(r, z) at d = AAclq = 2.93a (right) at the critical composition 
4> = 0.5, where there is no prewetting transition. Top: Those 
for (originally) weakly hydrophobic surface with Oq7 = 0.08 
and agao = 0.018 at x = 1.2 and x = 1-965. Here F int 
changes its sign on approaching the solvent criticality The 
surface becomes hydrophilic (</>' < 0) weakly for \ — 1.2 and 
strongly for x = 1-965 due to ionization. Bottom: Those 
for strongly hydrophobic surface with aj)7 = 0.2 and a^ao = 
0.01 at x = 1-2 and x — 1-965. Here the surface remains 
hydrophobic even for full ionization. 



3. Partial wetting on hydrophobic surface 

In the previous examples of hydrophobic colloids, the 
surface is in a completely wetted state above the precipi- 
tation line in the composition range 4> = 0.4 — 0.5. How- 
ever, for smaller <j>, a hydrophobic colloid surface can be 
partially wetted if 7 is not very large. 

In Fig. 16, at <j> — 0.3, we display such composition 
profiles (f>(r, z) around two hydrophobic colloids at d = 
44arj = 2.93a by changing x- The other parameters are 
common to those in Figs. 11 and 13. For X > Xp = 1-12, 
there appears a water-rich region partially wetting the 
colloid surface. The surface is partially wetted for \ = 
1.25, 1.55, and 1.85, but is completely covered by the 
water-like component at x — 2. Here, a in the non- 
wetted regions is 0.082, 0.058, and 0.051 for X = 1-85, 
1.55, and 1.25, respectively, while a is nearly equal to 
unity in the wetted surface regions. In these examples, 
the number of the counterions from the wetted surface 
is only about 0.4% of the number of the cations in the 
wetting water-rich region, where the latter is about 5600 
for x = 1.25. 

In Fig. 17, we vary d at <fi — 0.3 and x = 1-55 to see 



how the partial wetting of two hydrophobic colloids is 
changed between bridged and disconnected states. In 
the left panel, F; nt (d) is displayed as a function of d, 
where hysteresis is exhibited between these two states. 
In the right panel, we show <p(r, z) in the connected case 
at d = 60ao and in the disconnected case at at d — 72ciq. 
In the latter case, a water-rich droplet is partially wetting 
the left colloid. In these cases, a is 0.99 in the wetted 
part and about 0.05 in the nonwetted part. 



E. Antagonistic ions 

With antagonistic ions added, oscillatory response can 
arise against local disturbances even if the system is in 
a homogeneous state in the bulk region, as has been 
discussed in the linear theory in Sec. III. Here, we set 
9i = -92 = 13, a§7 = 0.2, a\a = 0.016, and 4> = 0.5. 

In Fig. 18, we show cj>(r) and U(r) as functions of r/ao 
around a single colloid, where three curves correspond 
to x = 1-961, 1.965, and 1.966. Oscillatory behavior is 
amplified with increasing x- The system tends to a ho- 
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FIG. 13: (Color online) <j>(r,z) (left) and U(r,z) (right) 
around two hydrophobic colloids above a precipitation line, 
where ajfr = 0.08, agoo = 0.018, and (j> = 0.425. Wetting 
layers of two colloids are bridged at d = 5.33a = 80ao (top) 
and are disconnected for 6.8a = 102ao (bottom) . The surface 
has become hydrophilic with <j>' < 0. 
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FIG. 14: (Color online) Normalized interaction free energy 
Fi n t(d)/T as a function of normalized distance d/ao above a 
precipitation line, where x = 1-94, 1.95 and 1.96 for three 
curves (from above). Wetting layers of two colloids are 
bridged for small d but are disconnected for large d. There is 
hysteretic behavior between these two states (see Fig. 15). 



mogeneous state far from the colloid for the smaller two 
X, but it is in a periodically modulated phase for the 
largest \, since the homogeneous state is linearly unsta- 
ble for x > 1-9651. In Fig. 19, we examine the case of two 
colloids, In its left panel, the interaction free energy Fj nt 
is plotted as a function of d for two colloids for \ = 1.95, 
1.961 and 1.965. In its right panel, the profile of <fi(r,z) 
is given for x = 1-965 at d = 49ao = 3.27a, where a 
homogeneous state is linearly stable. 
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FIG. 15: (Color online) Left: Normalized interaction free en- 
ergy Fi nt (d)/r, exhibiting hysteresis between bridged and dis- 
connected states of wetting layers for x — 1-95 in the box re- 
gion in Fig. 14. There can be a metastable state in an interval 
of d (left). Composition <j>[r,z) in wetting layers at the two 
states marked by the arrows in the left. 



Partial wetting above precipitation: <p — 0.3 d — 2.93a 




FIG. 16: (Color online) Partially or completely wetting pro- 
files of 4>(r, z) around two hydrophobic colloids for (j> — 0.3 
with d = 44a = 2.93a, where x = 1-25, 1.55, 1.85, and 2 
above the bulk precipitation value Xp = 1-12, The other pa- 
rameters are common to those in Figs. 11 and 13-15. 



VI. SUMMARY AND REMARKS 

In summary, we have examined how ionizable colloids 
influence the ion distributions and the composition field 
in binary polar solvents. These perturbations then give 
rise to the interaction free energy F[ nt (d) between two 
colloids as a function of the distance d between their 
centers. We summarize our main results, 
(i) In Sec. II, we have introduced a Ginzburg-Landau 
model in the presence of negatively ionizable colloids. 
The fundamental fluctuating variables are the composi- 
tion (/>, the ion densities n, , and the degree of ionization 
a, which are inseparably coupled in the presence of the 
selective solvation. Important parameters in the bulk 
free energy are the interaction parameter x (determined 
by the temperature T), the average composition </>, the 
average anion density no, and the solvation parameters 
cji. Those related to the colloid are the radius a, the 
molecular interaction parameter 7 representing the 
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Partial wetting above precpitation : 




FIG. 17: (Color online) Two partially wetted, hydrophobic 
colloids in a bridged or disconnected state for <j> — 0.3 at 
X = 1.55. The other parameters are common to those in 
Fig. 16. Left: Fint(d) /fcsT vs d/ao with hysteresis between 
the two states. Arrows indicate d/ao = 44, 60, 72. Right: 
4>{r,z) at d/ao = 60 (bridged) and 72 (disconnected). Profile 
at d/ao — 44 is in the top right panel in Fig. 16. 
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FIG. 19: (Color online) Left: Normalized interaction free en- 
ergy Fint/ksT as a function of d for two colloids with antag- 
onistic ions, where <2q7 = 0.2 > 0, a^ao = 0.016, and <f> — 0.5. 
Three curves correspond to x — 1-950, 1.961 and 1.965 from 
below. Instability occurs at x > 1.9651, so homogeneity is 
attained far from the colloids. Right: (j>(r, z) around two hy- 
drophobic colloids at d — 49ao = 3.27a for x = 1-965. 
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FIG. 18: (Color online) (j>(r) and U(r) around a hydropho- 
bic colloid for antagonistic salt with gi — ~g2 = 13, where 
a o7 — 0.2, a^ao = 0.016, and <j) = 0.5. Three curves cor- 
respond to x — 1-961 (green), 1.965 (blue) and 1.966 (red) 
Parameter M in Eq.(3.21) in the linear theory is equal to 
0.8808, 0.8344, 0.8224, respectively. Oscillatory behavior is 
amplified with increasing x- The system tends to be homo- 
geneous far from the colloid for the first two curves, while a 
mesophase is realized for the third curve. 



surface field, the density of the ionizable groups on the 
surface ctq, and the composition-dependent ionization 
free energy T(Ao — Ai</>), which appear in the boundary 
condition (2.24) on <fi and the mass action law (2.26) for 
a. 

(ii) In Sec. Ill, we have presented a linear theory of 
the electrostatic and compositional disturbances pro- 
duced by charged colloids, which is a generalization of 
the Debye-Hiickel and DLVO theories. In the linear 
scheme, the colloid interaction free energy Fi nt (d) 
is a linear combinations of two Yukawa functions 
e-<»( d - 2o )/(l + qi/a) 2 d as a function of the colloid 
separation distance d in terms of two characteristic wave 
numbers In the weak coupling linit g a — g\ — g^ — > 0, 
they tend to the DLVO interaction in Eq.(1.2) and the 
adsorption- induced attraction in Eq.(1.5). 



(iii) In Sec. IV, we have presented a theory of precip- 
itation on the colloid surface assuming a completely 
wetting layer, which is induced by the selective solvation 
of hydrophilic ions far from the solvent coexistence 
curve. This precipitation occurs at small x f° r relatively 
small 0~ (say, 0.1) (see the left panel of Fig. 4). However, 
the growth of the layer thickness is slow with increasing 
X for small ip. For <p > <p c precipitation occurs close to 
the solvent coexistence curve. 

(iv) In Sec.V, we have presented numerical results on 
precipitation and prewetting on the colloid surface for 
hydrophilic (7 < 0) and hydrophobic (7 > 0) surfaces 
in the nonlinear theory. We are particularly interested 
in the weakly hydrophobic surface satisfying Eq.(2.29). 
Such a surface is hydrophobic without ionization, but 
becomes hydrophilic with progress of ionization. Also 
the prewetting phase transition is more dramatic for 
such a hydrophobic surface than for a hydrophilic surface 
as in Figs. 6 and 8. Wetting layer formation occurs above 
a precipitation line, which weakly depends on the radius 
a and is located slightly above the bulk precipitation 
line x — Xp- Such layers undergo a bridging transition 
with a great change in the interaction free energy i*i nt 
as in Figs. 13-15. They either completely or partially 
wet the surface depending on the average composition <p 
as in Figs. 16-18. For antagonistic ion pairs, oscillation 
can be seen the composition and potential profiles as a 
function of the separation distance d as in Fig. 18, but it 
is largely masked in Fi nt in Fig. 19. 

We make some remarks. 
1) Our coarse-grained theory is inaccurate on the 
angstrom scale, but the solvation parameters gi and the 
ionization parameter Ai can be made very large, so the 
precipitation and prewetting transitions on the colloid 
surface have been predicted. We also note that the molec- 
ular volumes of the two components are often very dif- 
ferent in real mixtures. For example, those of D2O and 



18 



tri-methylpyridine (the inverse densities of the pure com- 
ponents) are 28 and 168 A 3 , respectively 37]. The coef- 
ficient C in Eq.(2.1) of the gradient free energy remains 
an arbitrary constant (39|, though we have set C = 2a^" 1 
inSec.V. 

2) We believe that the previous observations of colloid 
aggregation [6MTo[| should be induced by overlapping of 
enhanced adsorption or wetting layers on the colloid sur- 
face 14j. If the ions are neglected, colloid particles con- 
stitue a selective solute added in a binary mixture 0. 
Furthermore, if the selectivity is high, our previous the- 
ory [24[ indicates a solute-induced phase separation with 
a phase diagram as in Fig. 4. This aspect should be stud- 
ied in more detail. 

3) Wetting behavior remains unexplored in the presence 
of a highly selective solute. It becomes even more com- 
plex if the substrate itself is ionizable. We have realized 
both complete and partial wetting on ionizable colloid 
surfaces, but the information gained is still fragmentary 
because many parameters are involved in the problem. 

4) The effects of the critical fluctuations on the interac- 
tions between solid surfaces are very intriguing [29|-|32j. 
Ions should further promote bridging of highly adsorbing 
or wetting layers. 

5) For antagonistic ion pairs, the oscillatory behavior in 
the colloid interaction in Fig. 19 is rather mild, though 
it is evident in the composition and potential profiles. 
It becomes more evident in the interaction between two 
parallel plates, as in the case of liquid crystals (38j . 

6) For polyelectrolytes including ionized gels, there are a 
number of unsolved problems arising from selective sol- 
vation. Even in one-component solvents, ions interact 
differently with polymer segments and solvent molecules 
[23l | . In mixture solvents (water-alcohol), a wetting film 
should be formed around a chain, as stated in Sec. I 

23, 34 36] . The Manning- Oosawa counterion condensa- 
tion mechanism should be modified for mixture solvents. 
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Appendix A: Pair correlation functions 



1. Composition fluctuations 

We examine the structure factor of the composition 
fluctuations S(q) = {\(/>q\' 2 ) c , where <\>q is the Fourier com- 
ponent of 4>(r) with wave vector q and (■ ■ -) c denotes 
taking the thermal average. The mean-field structure 



factor reads [It 
1 



S(q) 



= T - T c 



1 - 



(Al) 



in terms of r in Eq.(3.6), r c in Eq.(3.9), and 7 P in 
Eq.(3.22). If the right hand side of Eq.(Al) is ex- 
panded with respect to q 2 , the coefficient in front of 
q 2 is C(l — 7p). Thus a Lifshitz point is realized at 
7 p = 1. (i) For 7 p < 1, S(q) is maximum at q = 
and S(q)^ 1 oc £~ 2 + (1 — 7p)g 2 for q <C k, so a thermody- 
namic instability occurs for r < r c at long wavelengths. 
To be precise, t c is the shift in this case, (ii) For 7 p > 1, 
S(q) has a peak at 

9m = (7 P - l) 1/2 «- (A2) 

The corresponding peak height is given by 

S(q m ) = C^n-'/lM 2 - ( 7p - l) 2 ]. (A3) 

which diverges as 7 P — 1 — > M = (k£) _1 . A mesophase (a 
charge-density-wave phase) should emerge for 7 p — 1 > 
M, as was observed experimentally |37| . 

Furthermore, from Eq.(Al), the quadratic equation 
(3.19) is identical to Siqy 1 = with q 2 = -k 2 A. In 
terms of q\ and (72 in Eq.(3.18), we obtain 



CS{q) 



(q 2 + ql)(q 2 + «§)' 



(A4) 



The inverse Fourier transformation of S(q) yields the 
pair correlation g(r) = {5(f>(r)S(j)(0)) o for the composi- 
tion fluctuations. It follows a sum of the two Yukawa 
functions, 



50) = 



(Ai - l)e-^ - (A 2 - l)e- 
4ttC(Ai - A 2 )r 



(A5) 



In particular, in the region \M — 1| < 7 P < M + 1, qi and 
qi are complex conjugates and g(r) behaves as 



9{r) 



-qnr 



4nCr 



cos(q/r) + (ql 



q] 



2 ^sin(g/r) 



) 



Zquqi 



(A6) 



where qn and qi are are given in Eq.(3.24). 



2. Ion fluctuations 

We eliminate the composition fluctuations assuming 
their Gaussian distribution, where the ion densities are 
held fixed. The resultant ion-ion potentials read [3, E3] 



Vij(r) = ZiZj 



T9i9j 1 r/i 



er AttC r 



(A7) 



where Z\ and Z 2 are ±1 in the monovalent case. The sec- 
ond term is the composition-induced interaction decaying 
exponentially with the correlation length £ = (C/r) 1 / 2 . 
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It is attractive among the ions of the same species (i = j). 
It dominates over the Coulomb repulsion for g 2 > AttC£b 
in the range r < £, under which there should be a 
tendency of ion aggregation. In the antagonistic case 
(<7i32 < 0), the cations and anions repel one another 
for |gi<?2| > 4:TtC£b in the range r < £, leading to charge- 
density- wave formation near the criticality. Note that the 
shifted correlation length £ in Eq.(3.11) has appeared in 
the colloid-colloid interaction, where both the composi- 
tion and ion densities are eliminated. 

In our recent review papers [T§l . [l9j , we have further- 
more calculated the structure factors among the ion den- 
sities Gij(q) = (niqn'q) Jno- Further using Eq.(A4) we 
find 



Gu(?) = l- 



1/2 



7 2 [m + w(u + 1)] 2 /2 



tt+1 (u+l)(u + X 1 )(u + A 2 ) : 



Gi 2 (g) 



1/2 7 2 [u; 2 (u + l) 2 -w 2 ]/2 



+ 



u + 1 (u + l)(u + Ai)(u + A 2 )' 



(A8) 



where u = q 2 / 'k 2 and w = (gi + g-2)/{gi — .92)- The 
G 22 {q) is obtained if w in Gu(q) is replaced — w. The 
inverse Fourier transformation of these structure fac- 
tors gives rise to the pair correlation functions gij(r) — 
{5ni(r)5nj(0)) c /n . We notice that the terms pro- 
portional to e~ KV /r cancel to vanish from the relation 
(1-Ai)(l-A 2 ) =7p- Thus, 

9ij (r) = S(r)S l3 - (Kfe~^ + K$ e -* r )±, (A9) 

where the 5 function appears due to the self correlation 
and K-y and K 2 are appropriately defined constants. 

Appendix B: Calculations in the linear theory 
for one and two colloid particles 

First, we seek the fundamental profiles Ua(r) and tpo( r ) 
around a single colloid particle induced by by the bound- 
ary conditions (3.12) and 83.13), where the surface charge 
is fixed. They depend only on the distance r from the 
colloid center. For r > a we may set 



Uo(r) =g a (a 1 e~^ 



a 2 e 



-qir 



)/r, 



Mr) = {Pie- qir -^- q2r )/r- 



(Bl) 
(B2) 



For r < a, we have Uo(r) = Uq(o) =const. from V 2 /7o = 
0. From Eqs.(3.7), (3.8), and (3.21), we obtain 

ft = (1 - Ai)oi, (B3) 
[A,-M 2 + 7 2 ]ft = 7 2 a,, (B4) 

which hold for i = 1,2. The boundary conditions (3.12) 
and (3.13) give 

a[ -a' 2 = -MBa/ga, (B5) 

4(l-Ai)-a 2 (l-A 2 ) = -7/C, (B6) 

where a'- = a^(l + <^a)e~ 9ia /a 2 . Using the relation (Ai — 
1)(A 2 — 1) = 7p, we solve these equations to obtain 

e* a Bi(Ai-A 2 ) 



(l + gjoXl-Ai)' 



(B7) 



where B\ and B 2 are defined in Eq.(3.26). We thus con- 
firm Eqs.(3.25) and (3.26). In this one-colloid case Af2 
in Eq.(3.17) is written as -F se if- Some calculations give 



2ira 6 T 



1 + q\a 1 + q 2 a 



(B8) 



where E\ and E 2 are defined in Eq.(3.31). 

Next we consider two colloid particles at positions v\ 
and r 2 separated by d = \r± — r 2 | > 2a under the con- 
dition of fixed surface charge. In the colloid exterior 
(|r — ri| > a and \r — r 2 \ > a), U and ijj are expressed as 

U = U (\r - n|) + E7b(|r - r 2 \) + W(r), (B9) 
%l> = M\r - ri\) + M\r - r 2 \) + 9{r), (BIO) 

where Uq and tpo are the fundamental profiles for a single 
colloid. In the colloid exterior, we expand the corrections 
W and "J around the center of the first colloid at r\ as 

W = g a J2 a Mn\r - n\)P t (0i), (Bii) 
* = 5^(1 - Xi)aak t { qi \r - ri\)Pi{B x ), (B12) 

i,t 

where an [i = 1,2 and £ — Q, 1, - • - ) are unknown co- 
efficients to be determined below. The Pi(9\) are the 
spherical harmonic functions with 9\ being the angle be- 
tween r — T"i and r 2 — r\. We introduce the modified 
spherical Bessel functions it(x) and kt(x) [49]. They sat- 
isfy i'J + 2i' e /x-£(e+l)ie/x 2 = u and k' e ' + 2k' e J x - £(£ + 
l)ki/x 2 = ki, where we write i\ = di^/dx, k' f — dkg/dx, 
i" t = d 2 ii/dx 2 , and fc" = d 2 k(/dx 2 . We have u ~ x e as 
x — > and kg ~ e~ x as x — > 00. In particular [50], 

zo(x) = sinhx/a;, feo(aO = e~ x /x. (B13) 

Thus, i^Pf and /c^-Pf satisfy the Hehnholtz equations, 

(V 2 -qf)[i i (q i \r-r 1 \)P i (e 1 )]=0, 

(V 2 - qf)[k e ( qi \r - ri\)P t (6i)] = 0. (B14) 

With these relations and Eq.(3.7), we derive Eq.(B12) 
from Eq.(Bll). In Eqs.(B9) and (BIO) we also need to 
expand Uq(\v — r 2 \) and 4'o{\ r ~ r 2\) around r\ in terms 
of Pi(8i). To this end, we use the following mathematical 
relation [SEH, 



,-g\r-r 2 \ 



r - r 2 \ 



Y^^+^QhiqdMqlr-r^P^), (B15) 



which holds for Re(g) > and in the region \r — r%\ < 
d = ri — r 2 \. On the other hand, in the interior of the 
first colloid \r — r\\ < a, we have V 2 U — so that we 
may assume the expansion, 



U = U Q (a)+g a Y,bi\r-r 1 \ e P e 



(B16) 
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We can calculate the coefficients an and bi from the 
boundary conditions (3.12) and (3.13) and the continuity 
of U at the colloid surface. We are interested in the free 
energy deviation Afi in Eq.(3.17). For two symmetric 
colloids, we obtain 



we obtain simple Yukawa forms, 



AQ/T = 47ra 2 [a(V>)o-7(£/}o], 



(B17) 



where (• • -)o denotes taking the surface average on the 
colloid 1 (= f* d0i sin 0i (■ • • )/2 at \r - n| = a). Thus 
Afl arises from the terms with I = in Eqs.(Bll), (B12), 
and (B15). For I = 0, the boundary conditions (3.12) and 
(3.13) simply yield 

^2 <li[ a ioQiko{qid)i' i0 + a i0 k' i0 ] = 0, (B18) 

i 

^(1 - \i)qi[oti qiko(qid)i' i0 + a i0 k' i0 ] = 0. (B19) 

i 

For simplicity, we write io(qi(i), i' (qia), ko(qia), and 
k' (qia) as i^, i' i0 , kio, and k' i0 , respectively, suppress- 
ing qia. In Eq.(B18) there is no contribution from the 
electric field within the colloid (oc e p ). This is because 
the angle average of (r — r\)-E vanishes from Eq.(B16). 
For each i, it follows the relation, 

a i0 = -a l q i k {q l d)i / l0 /k' l0 . (B20) 
Elimination of yields 
(U)o = U (a) + g a ^2uiko(qid)[i l0 - k M i' i0 /k' i0 \. (B21) 

i 

The expression for (tp) also follows in the same manner. 
Further, using the relation, 



i (x) - i' (x)ko{x)/k' (x) = e x /(l + x), 



(B22) 



-qi (d—a) 



(U) = U (a) +9a j:a iJ — j - d , (B23) 

U(d- 
ljd~ 



(~\ _ \.\ p —Qi(d—a) 

Wo = *>(«) + !> {1 + m y ■ (B24 > 



We also have bo = [(U)o — Uo(a)]/g a from the continuity 
of U. Substitution of Eqs.(B23) and (B24) into Eq.(B17) 
leads to the interaction free energy 



F ilrt = Af2 — 2R 



self J 



(B25) 



given in Eq.(3.30). 

Finally, we calculate the terms with t > 1, though 
they do not contribute to Ail in the linear theory. From 
Eqs.(3.12) and (3.13) we express an in terms of bt and 
Qj as 

a ie = [eb e a e ~ 1 ^ l e p /e - arffkeiqid^/qik'n, (B26) 
where = (1 - A 2 )/(Ai - A 2 ) and fi 2 = (1-Ai)/(A 2 -Ai) 
in the first term. We write in — ie(qia), — i'^qia), 
ku — ki(qia), and k' if = fc^(g,a). Requiring the continuity 
of the potential, we determine bi in the form, 



ha 



e-i _ Yh a-iqikt{qid)r]u 

a - e(e p /e)J2iHiWn/qi 



(B27) 



where r)u = %u - i' u k it /k' u and 



knjk' it The term 



proportional to e p /e in the denominator in the right hand 
side arises from the boundary condition (3.12). 
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